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^ I We present precision calculations of the processes e+e^ -^ 4 fermions, in which the double 

■^ ■ resonant W^W~ intermediate state occurs. Referring to this latter intermediate state 

p ■ as the signal process, we show that, by using the YFS Monte Carlo event generators 

ji. . YFSWW3-1.14 and KORALW 1.42 in an appropriate combination, we achieve a physical 
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The award of the 1999 Nobel Prize for physics to G. 't Hooft and M. Veltman, and the 
success of the predictions of their formulation [0] of the renormalized non-Abelian quantum 
loop corrections for the Standard Model of the electroweak interactions in confrontation 
with data of LEP experiments, underscores the need to continue to test this theory at 
the quantum loop level in the gauge boson sector itself. This emphasizes the importance 
of the on-going precision studies of the processes e^e~ — >• W^W~ + n(7) — >• 4/ + 71(7) 
at LEP2 energies 0-|^, as well as the importance of the planned future higher energy 
studies of such processes in LC physics programs j^-^. We need to stress that hadron 
colliders also have considerable reach into this physics and we hope to come back to their 



roles elsewhere 10 



In what follows, we present precision predictions for the event selections (ES) of the 
LEP2 MC Workshop lll|, for the processes e+e~ -^ W^W~ + n{^) ^ 4/ + n(7), based 



on our new exact 0{a)prod YFS-exponentiated LL 0{a ) FSR leading-pole approxima- 



tion (LPA) formulation, as it is realized in the MC program YFSWW3-1.14 |[T^,|T^, in 
combination with all four-fermion processes MC event generator KoralW-1.42 [H, T^] so 
that the respective four-fermion background processes are taken into account in a gauge- 
invariant way. Indeed, gauge invariance is a crucial aspect of our work and we stress that 
we maintain it throughout our calculations. Here, FSR denotes final-state radiation and 
LL denotes leading-log as usual. 

Recently, the authors in Refs. [1^] have also presented MC program results for the 
processes e^e~ -^ W^W~ + ^(7) -^ Af + ^(7), n = 0,1, in combination with the 
complete background processes that feature the exact LPA 0{a) correction. Thus, we 
will compare our results, where possible, with those in Refs. fl^ in an effort to check the 
over-all precision of our work. As we argue below, the two sets of results should agree at 
a level below 0.5% on observables such as the total cross section. 

More specifically, in YFSWW3-1.13 |T^, the leading-pole approximation (LPA) is 



used to develop a fully gauge- invariant YFS-exponentiated calculation of the signal process 
e~^e~ — * W'^W~'+n{'~f) -^ Af + n{'j), which features the exact 0{a) electroweak correction 
to the production process and the C(a^) LL corrections to the final-state decay processes. 



The issue is how to combine this calculation with that of KoralW-1.42 in Ref. [Il4| , p!5[| for 
the corresponding complete Born-level cross section with YFS-exponentiated initial-state 
0{a^) LL corrections. In this connection, we point out that the LPA enjoys some freedom 
in its actual realization, just as does the LL approximation in the precise definition of the 
big log L, without spoiling its gauge invariance. This can already be seen from the book 
of Eden et al. |]T^ , wherein it is stressed that the analyticity of the S-matrix applies to the 
scalar form factors themselves in an invariant Feynman amplitude, without any reference 
to the respective external wave functions and kinematical (spinor) covariants. The classic 
example illustrated in Ref. [llTl is that of pion-nucleon scattering, with the amplitude 



M = u{p2)[A{s,t) + B{s,t){4,+ 42)]u{pi), (1) 

where the Pi are the nucleon 4- momenta, the g, are the pion 4-momenta, u{p) is the usual 
Dirac wave function of the nucleon, and the invariant scalar functions A{s,t) and B{s,t) 
of the Mandelstam invariants s = {pi + qi)"^ , t = {q2 — qiY realize the analytic properties 



of the S-matrix themselves in the complex s and t planes. This means that, whenever we 
have spinning particles, we may focus on the analogs of A and B in eq. (|l]) in isolating 
the respective analytic properties of the corresponding S-matrix elements. We note that 
Stuart |T^ has emphasized this point in connection with the production and decay of 
Z-pairs in e~^e~ annihilation and in connection with the production and decay of single 
VT's in e"'"e~ annihilation. What this means is that, in formulating the Laurent expansion 
of the S-matrix about its poles to isolate the dominant leading-pole term (the LPA is then 
realized by dropping all but this leading term), we may focus on only A and B, or we 
may insist that in evaluating the residues of the poles in the S-matrix the wave functions 
and kinematical covariants are also evaluated at the pole positions. When we focus only 
on the analogs of A and B in formulating the LPA, we shall refer to the result simply as 
the LPAq; when we also evaluate the wave functions and kinematical covariants at the 
pole positions in isolating the poles in the analogs of A and B for the LPA, we shall refer 
to the respective result as the LPA^. As Stuart stressed as well, both the LPA^ and the 
LPAb are fully gauge-invariant. 



For the process under discussion, a general representation is [1 



M = J2^j^ji{lkQi}), (2) 

j 

where {ij} are a complete set of kinematical covariants which carry the same trans- 
formation properties as does Ad, and the Lorentz scalars {qkQi} are a complete set of 
Lorentz scalar invariants for the external 4- momenta oi Ai. In the LPA^, we make 
a Laurent expansion of the Aj and retain only their leading poles, without touching 
the {ij}] in the LPA?,, we also evaluate the ij at the position of the respective lead- 
ing poles. Evidently, in the latter case, we must make an analytic continuation of the 
phase-space point originally associated with the {£j} to a corresponding such point for 
the respective pole positions. See Ref. [1^ for an illustration of such a continuation in 
the context of the YFS-exponentiated exact 0{a) calculation for the production process 
in e~^e~ — > W^W~ +n{'y) -^ 4/ + n{'y), and Refs. for a similar illustration in the con- 
text of the 0{a) correction to e'^e' -^ W^W~ — * 4/. Having isolated the appropriate 
realization of the LPA at the level of A^, it must still be decided whether to treat the 
phase space used to integrate the cross section exactly or approximately to match what 
was done for the {ij} in the case of the LPA;,. In all of our work, we stress that we always 
treat the exact phase space, both in the LPA^ and in the LPA^. 

In the context of YFS exponentiation, we realize the LPA as follows, as was briefly 
described already in Ref. |1^. Taking the respective 4-fermion plus n-photon process 



kinematics to be as given by (here, (ir„+4 is the respective phase space differential with 



the appropriate normalization): 

e-(pi) + e+fe) - fiin) + /2(r2) + f[{r[) + f^{r',) + 7(^1), ...,7{K) 

<^n = -Ti / dTn+4{pi + P2] n, Ts, r[, r'^, h, ..., kn) 

flux J 

Yl 5Z \^tfiPi^P2,ri,r2,r[,r2,ki,...,kr, 



(3) 



r, 



Ferm. Spin Phot. Spin 

and that of the corresponding W^W~ production and decay process to be as given by 

e-(pi) + e+(p2) ^ W-{q{) + W+iq^), 
W-{q,) ^ /i(n) + /2(r2), ly+fe) - /{ W) + /^(r^), 
1 



cr„ 



rfrn+4(Pi + P2; ri, r2, r^, r2, fci, ..., /c„) (4) 



5Z 5Z l-^ipA(Pi>P2,ri,r2,r'i,r2,/i;i,...,/c„)|2 

Ferm. Spin Phot. Spin 

in the context of YFS exponentiation [|19|,|^, we proceed according to Refs. [p!^ ,p!9|,pO| 



Mff{pi,P2,ri,r2,r[,r2,ki,...,kn) M^^l^{pi,p2,ri,r2,r[,r2,ki, ...,kn) 

^ X] ^Prod^\Pl^P2,(ll,(l2,kl,...,ka) 



Phot. Partitions 



X -^r— M%^^ly^^{qi,ri,r2,ka+l, ...,h) X TT^^ -MoeealAafe, r^, r^, Vi, ..., fe 



D(g,) = g,2 - M\ M^ = (M^ - tTwMw){l - T^/M^ + 0(a')), 
gi = ri + r2 + fca+i + ... + h; g2 = ^'1 + ^2 + Wi + ••• + ^n, 

so that M^ is the pole in the complex q^ plane when q is the respective W 4-momentum, 
and My/ and Tw are the on-shell scheme mass and width, respectively. The residues in 
(^ are all defined at q]^ = M^ with a prescription according to whether we have LPA^ or 
LPAfo, so that d^) is our YFS generalization of the formula in eq. (12) in the first paper 
in Ref. 0: 

, , ,,-^ A+ (Ml) A7 (M2) 

A<(")=5^n,„,,(Mi,M2)^^^^, n = 0,l, (6) 

where Di = D{qj) and Mf = M^. We stress that, unlike what is true of the formula in 
eq. (12) in the first paper in Ref. [^ and in eq. (P) here, in eq. (|^) n is arbitrary. The 
sum over "photon partitions" is over all 10*^ possible attachments of n photons to the 
six external fermion lines and the two W"^ lines (one for the W production and one for 
the W decay, respectively). We make the further approximation that M^ = M^ in the 
residues in (|^), always maintaining gauge invariance, as explained. Equations (3) and (4) 



in Ref. |2^ then give us, in the presence of renormahzation-group-improved perturbation 
theory, for the representation 



•^ip^bi'^s, ri, r2, r[, r'^, h, ..., K) = J^ M'f\puP2, n, rs, r[, r'^, h, ..., A;„), (7) 

i=o 

where Ad"" is the j-th virtual photon loop contribution to the residues in Ai^py^, the 
identifications 



,(n) 



^ (Pi,P2, n, r2, r'l, r2, fci, ..., kn) 



r=0 



[aB'Y 



Tl 



(8) 



where 5' is now, for the LPA^ case to be definite, the on-shell virtual YFS infrared 
function, which reduces to that given in eqs. (8) and (9) in Ref. [^ when we restrict our 
attention to the production process, and a is indeed a(0) when it multiplies B' here. Let 
us keep this limit of B' in mind, as we focus on the gauge invariance of YFSWW3-l.il in 
Ref. [|12[, which treats the radiation in the production process, and on that of YFSWW3- 
1.13 and YFSWW3-1.14, in which the radiation from the decay processes is also treated. 
In the LPAa case, the corresponding B' function is off-shell. Let us discuss first the LPA^ 
case and comment later on how the corresponding results for the LPA^ case are obtained. 
Here, since the SU{2)l x f/(l) Ward-Takahashi identities require (see eq. (47) in 
Ref. m 



k^'Ml = 0, 



k^'Mt 



T^M^, k^M"^^ = ±./JI^M'^' 



(9) 



for fiv denoting the squared V boson mass (so that, for V = W, nw = ^^), we find 
that B' is SU{2)l x [/(l)-invariant from the equations in (|D and our result eq. (8) in 
Ref. ||20|. From eq. (|^) it then follows that the infrared residuals m 
f/(l)-invariant. Here, x ^"^^ 
renormalizable gauges and we use the notation of Ref. ^^ , so that M^ is their respective 
amplitude for the emission of a Z of Lorentz index /i and 4-momentum fc, and M'^ is their 
corresponding amplitude for the emission of a x with the same 4-momentum, etc. 
Introducing eq. (|) into (H) gives 



'.^•_^ are also SU{2)l x 
^ are the usual unphysical Higgs fields in our general 



■^tpA,iP^^P2, n, r2, r[, r2, /ci, ..., kn) = e"^' ^ m'f\puP2, ri, r2, r[, r'^, /ci, ..., kn). (10) 



Equation (2.13) in Ref. [T^ and eq. (7) in Ref. [^ then give our n-photon differential 
cross section, for P = pi + P2, P = 0, as 



d'^LPAb — ^ 









(n) 



n'=0 



d^ri d?r2 d?r[ d^r'2 
'1 '2 "1 "2 



'111 



where we note that, when we only focus on the production process in eq. (0), R is the 
produced WW intermediate state; R = ri + r2 + r[ + r'2. Using the second theorem of the 
YFS program (eq. (2.15) in [jl9[), "we get 



(n) 



E"^-' 



n'=0 



S{h) ■ ■ ■ S{k^)Po + Yl "^(^l) ■ ■ • '^(^.-l)'^(^.+l) ■ ■ ■ S{kn)Pl{k, 



i=l 



(12) 



+ ■ ■ ■ + y ^ j[kijPn-l{ki, ■ ■ ■ , ki-i, /Cj+i, ■ ■ ■ , kn) + Pn[kl, ■ ■ ■ , kn), 



i=l 



where the real emission function S{k) is given by Sprod{k), the real emission infrared 



function in eq. (8) in Ref. [2C] for on-shell W's, when we only focus on the emission 

(13) 



from the production process as we did in Refs. |T^,|13[. Since, in general 

S{k) = Sprod + Soeci + SDec2 + Sint, 

with 



Sprod{k) 



a 
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(14) 
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(15) 



Suec^ik) 



a 



4772 



Qmh-^] +Q'^Qw(h 



Q'2Qw I -rr 



kr'^ 
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M2 Y 

kAq2j 
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Slnt{k) 



Att"^ V kpi kp2 kri kr2 



+ «d7-<3.i4 
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S 



Prod 'JDeci ^Dec2 



Sr 



(16) 



(17) 



is really composed of the scalar product of emission currents {Jh{k)} with k^jj^{k) = 0, 



eq. 



]) is also SU{2)l x [/(l)-invariant. Here 

Aqi = analytical continuation of g^ to the point qf 



Ml- 



'li 



This analytical continuation, already described in Ref. [|1^], does not spoil the gauge 
invariance, as we see from Eqs. (|13|,|18D. It follows that the hard-photon residuals {/5„} 
are also SU{2)l x f/(l)-invariant. 

Substitutung ([T^) into (0) we finally get the SU{2)l x f/(l)-invariant expression, 
which is the fundamental formula of our calculation, 



da 



LPAt 



23iaB'+2aB 



1 



{Any 



^4 i|/{pi+P2-gi-g2)+ 



D 



n ,sr^ Q i^j -iyk 

n=l '^i 



^ Pn[ki, ..., kn) 



(fri d?r2 d^r[ d^r!^ 
^1 ^2 r'-i rL 



(19) 



where we have defined 



D 



—S 
ko 



'iy-k 



9{Krna. ' \k\) 



2aB 



d^k 

kn 



0{Kma. - \k\)~S{k). (20) 



This shows that the parameter K^ax '^ \fs is a dummy parameter on which da does 
not depend. In (|1^), when the complete value of S{k) is used, then all W^ radiative 
effects are contained in the respective /3„ residuals, in accordance with the YFS theory 
in Ref. [|1^, as the non-zero widths of the W^s prevent any IR singularities when a W 
radiates a photon in @). In our work in YFSWW3, as we indicate below, we make the 
approximation of dropping all interference effects between the production and decay stages 
and between the two decay stages of (H). This means that we drop the Sjnt{k) in S{k) 
in (p!3D and in (|19|), so that the YFS theory then determines the corresponding forms of 
the YFS functions /5„, B' and D as also having the respective interferences dropped. This 
approximation, which resums a certain class of large W radiative effects, corresponds to 
the YFS exponentiation of the W production and decay radiation in the LPA, neglecting 
of all interferences between the production and decay stages and between the two decay 
processes. 

We will now comment further on our use of ([l9|). Since the residuals on the RHS 
(right-hand side) of (|^) are on-shell amplitudes in the SU{2)l x f/(l) theory, for both the 
production and the decay process, it follows that they satisfy the renormalization group 
equations |^,0 for the SU{2)l x f/(l) theory, as explained in Ref. [pl|]: 



^| + /?.(te})^ 



leM9iR})m^ 



d 



jR 



dnij 



R 



iviigrR}) r = 



(21) 



for r = A^ 



(n)Ai,A2 
Prod 



,M 



in') 



where fi is the arbitrary renormalization point, {gin} are 



the respective renormalized SU(2)i x U{1) couplings, and {rrijij} are the corresponding 



renormalized mass parameters, etc., as defined in Ref. |2^. It follows tliat any two 
scliemes for computing M-pl^d ^i-^DecX- ^^^ related by a finite renormalization-group 
transformation. Tlius, tlie complex pole scheme (CPS) and the fermion loop scheme 
(FLS) [^ values oi M.pl^^' ^ i-^Dec x ^^^ related by such a finite renormalization-group 
transformation. Specifically, if we denote CPS and FLS values of any quantity A by 
A(CPS) and A(FLS), respectively, then we have the identity 

where Tun is the respective unrenormalized value of F and the Zr(R), R = CPS, FLS are 
the respective field renormalization constants. It therefore follows that we have the finite 
renormlization group transformation 

r(CPS) = Zr(cPS)^f(FLS)r(FLS) (23) 

connecting the FLS and CPS schemes. Of course, as the usual implementation of the FLS 
scheme omits a gauge invariant set of contributions to the heavy vector boson widths, for 
example, in checking the the result in (^) we must omit the corresponding contributions 
on both sides of the equation for consistency, as is usualy the case when we compare 
renormalzation group improved quantities. Indeed, as we take the normalization points 
for the two schemes to be the complex pole position M^ = ^w, and as we evaluate 
■^Prod ^1-^Dec A ^^ ^^^^ normalization point, the only difference in using the FLS in- 
stead of the CPS will be the approximate treatment in the FLS of the actual values of 
■^Prod ^ i-^Dec A ' ^^"^ /^VK fo^ the polc positiou, as already shown in eq. (8) of the first 
paper in Ref. [Q: thus, for example, if we keep only fermion loops in 0(a), only the 
lowest-order width appears in eq. (6) of the first paper in Ref. ; thus if one would ex- 
press the resulting ^w = M^ in terms of the respective on-shell mass and width, only the 
lowest-order part of Tw would be given properly. Similarly, as both schemes normalize at 
/iiy, the difference in the residues Aipl^^' ^ -i-^oec a i^ t^sX in the FLS only the fermion- 
loop contributions are retained, whereas the CPS keeps all loops. Evidently, we may 
extend our 0{a) calculation in the FLS by using the complete value oi nw and including 
all the one-loop corrections and attendant 0{a) real corrections from Ref. [PB|, as we did 



in Ref. [|T2|. We conclude that, after adding in the entire 0{a) correction from Ref. p6| 



our LPA exact 0{a)prod YFS-exponentiated calculation arrives at the same amplitudes, 
independent of whether we started with the FLS or the CPS. 

All of the above results extend directly to the calculation when we use LPAq ampli- 
tudes, as these are also gauge-invariant by the gauge-invariance of our leading poles in 
the S-matrix. Thus, the only change we must make is that the respective residues must 
be calculated in the LPA^ rather than in the LPA^; for example, in eqs. ([T^^^, g^ would 
be used instead of Aqi, etc. We have done this, as we further illustrate in the following. 

Let us now comment on the issue of the pure FSR YFS exponentiation for the decay 
processes treated in the LPA. We proceed in analogy to what is done in Ref. [^ for 



the MC YFS3 for the respective FSR. Specifically, for both decay residue amplitudes 
•^Dec A ' ^^ ^^y have contributions to the respective hard-photon residuals /3„ due to 



emission for the final-state decay processes; in these we follow the procedure, described in 
Refs. P3,p8l and already illustrated in Ref. p9[, for including these contributions, using 



the same YFS methods as we used above for radiative effects to the initial, intermediate 
and final states. Here, we shall neglect all interference effects between the production 
and decay processes as we explained above; this is analogous neglecting all interference 
effects between the initial and final states in the f3n in Refs. p7|,|29[]. We note that it 
is possible to retain these interference effects, as we have illustrated in the exact 0{a) 
YFS-exponentiated BHWIDE MC in Ref. [^ for wide-angle Bhabha scattering and, more 



recently, using the new CEEX exponentiation theory in Ref. [^ , to all orders in a in the 
new KK MC |2^,^ for the 2 fermion processes from the r threshold to 1 TeV. In this 
way we see that the use of eq. (|19D to include exponentiation of the FSR is fully realizable 
by Monte Carlo methods we already tested. We stress that, for the same reasons as we 
gave for the exponentiation of the complete process, these FSR contributions to the /?„ 
are fully gauge-invariant. 

In the current version of YFSWW3, version 1.14, we also drop the SDea terms in S 
in eq. ([T3|) and the corresponding terms in the functions /?„, B' and D, and include FSR 
using the program PHOTOS [0, which gives us a LL O(a^) realization of the FSR in 
which finite pt effects are represented as they are in the 0{a) soft-photon limit. This LL 
implementation of FSR is fully gauge- invariant. The ratio of BRs is then used to obtain 
the 0{a) correction in the normalization associated with the 0{a) correction to the de- 
cay processes themselves. Evidently, these ratios of BRs are also gauge- invariant. As we 
illustrate below, for the corresponding non-factorizable corrections we use an efficient ap- 
proximation in terms of the so-called screened Coulomb ansatz [^], which has been shown 



m 



m 



to be in good agreement with the exact calculations for singly inclusive distributions p3 
This ansatz is gauge-invariant. 

We also point out that the current version 1.14 differs from version 1.13 in Ref 
that it uses a different renormalization scheme. Specifically, the scheme used in version 
1.13 is the so-called G^ of Ref. ||26|, in which the weak-scale coupling p6[ acp is used for 
all terms in the virtual correction, except those that are infrared-singular, which are given 
the coupling a = a{0). In the renormalization-group-improved YFS theory, as formulated 
in Ref. [^, all the terms in the amplitude that involve corrections, in which the emitted 
photon of 4-momentum k has fc^ —>■ 0, should have the coupling strength corresponding to 
a{0) - not just those that are IR singular. We therefore have introduced into YFSWW3 
this requirement of the renormalization-group-improved YFS theory to arrive at version 
1.14. We refer to this scheme as our scheme [A). According to the renormalization-group- 
improved YFS theory, it gives a better representation of the higher orders effects than 
does the G^ scheme of Ref. [^. We stress that this scheme (A) is also gauge- invariant. 
The main effect of this change in renormalization scheme between versions 1.13 and 1.14 
is to change the normalization of version 1.14 by ~ —0.3% to —0.4% with respect to that 
of version 1.13 



11 



The generic size of the resulting shift in the YFSWW3 prediction, which we just 
quoted, can be understood by isolating the well-known soft plus virtual LL ISR correction 



to the process at hand, that has in 0{a) the expression [^] 

SltkLL = /3 In fco + ^ Ql + ^ - 2) , (24) 

where P = —{L — 1), L = ln(s/r7ig), and k^ is a dummy soft cut-off that cancels out of the 
cross section as usual. In the G^ scheme of Refs. [^, which is used in YFSWW3-1.13, 



only the part P \n kQ + {a / it) {n'^ / 3) of (5^^^^^ has the coupling a(0) and the remaining part 
^^^IsRLL ^^^ ^^^ coupling ac^ = a;(0)/(l — 0.0371). The renormalization-group-improved 
YFS theory implies, however, that a(0) should be used for all the terms in S^^^ j^^. This is 
done in YFSWW3-1.14 and results in the normalization shift ((q;(0) — aGf,)/Ti) (1.5L — 2), 
whichis ~ —0.33% at 200 GeV. This explains most of the change in the normalization of 
YFSWW3-1.14 with respect to YFSWW3-1.13. Moreover, it does not contradict the 
expected total precision tag of either version of YFSWW3 at their respective stages of 
testing. We stress that, according to the renormalization-group equation, version 1.14 
is an improvement over version 1.13 in that it better represents the true effect of the 
respective higher-order corrections. 



For the purposes of cross checking with ourselves and with Ref. [16|, we also created a 



second scheme, scheme (B), for the realization of the renormalization in YFSWW3-1.14. 



In this scheme, we put the entire 0{a) correction from Refs. ||26| at the coupling strength 
a = a{0). Since the pure NL hard 0{a) correction is only ~ —0.006 at 200 GeV, scheme 
{B) differs in the normalization from scheme [A) by ~ {a{0)/aG^ — 1)(— 0.006) = 0.0002, 
which is well below the 0.5% precision tag regime of interest for LEP2. Thus, scheme 
(B), which is used in Ref. [^, is a perfectly acceptable scheme for LEP2 applications. It 
gives us a useful reference point from which to interpret our comparison with the results 
of RacoonWW from Ref. [1^, which we discuss below. 



Having presented our gauge-invariant calculation as it is realized in the MC YFSWW3- 
1.14, we will now turn to illustrating it in the context of LEP2 applications. Specifically, 
we always have in mind that one will combine the cross section from YFSWW3 with that 
from KoralW-1.42 |T3,|T5[ MC, which is capable to calculate the non-resonant background 



contribution in a gauge-invariant way. We can do this in two ways, which we will now 



briefly describe and refer the reader to Refs. ^^ for the more detailed discussion. In 
the first way, we start with LPAq and we denote the corresponding cross section from 
YFSWW3-1.14 as cr[Ya). It is corrected for the missing background contribution by 
adding to it a correction Aa{K) from KoralW-1.42 to form 

aY/K = cTiYa) + AaiK), (25) 



where Aa{K) is defined by 



Aa{K) = a{K^)-a{K,). (26) 



Here, the cross section o"(i^i) is the complete 4-fermion result from KoralW-1.42 with 
all background diagrams and with the YFS-exponentiated 0{a^) LL ISR and the cross 
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section a{K^) is the restricted CC03 Born-level result from KoralW-1.42 - again with the 
YFS-exponentiated 0{a^) LL ISR. The result in eq. (H) is then accurate to C(f ^). 
Alternatively, one may start with the cross section for LPA^ ;, in YFSWW3-1.14, which we 
refer to as cr(Ya) and cr(Yb) correspondingly, and isolate the respective YFS-exponentiated 
0{a) correction, Acr(y), which is missing from the cross section cr(-ft'i) as 

Aa{Yj)=a{Y^)-a{Y,), (27) 

where criY^) is the corresponding cross section from YFSWW3-1.14, with the non-leading 
(NL) non-ISR 0{a) corrections to /?„, n = 0, 1, switched off. Then the cross section 

aK/Y = <y{Ki) + Aa{Yj) (28) 

has the accuracy of 0{^jj^). We have checked that the results ay/K and ctk/y are 
numerically equivalent at the 0.1% level of interest to us here. In the following we only 
show results from the former. For completeness, we also note that we sometimes identify 
cr(Yi) = CT{Ya)^ '^(^2) = o"(^b), cr(^3) = (^(-^s)? and a{K2) is to be identified as the cross 
section from KoralW-1.42 with the restricted on-pole CC03 Born-level matrix element 
with YFS-exponentiated 0{a'^) LL ISR. This latter cross section is a future option of 
KoralW |jlO|. It would allow further combinations of YFSWW3 and KoralW with the 
desired 0{^j^) accuracy. Such combinations would be of use in cross checks of our 
work. 

We now illustrate our precision predictions using ay/K- We have checked that the 
correction Aa{K) is small, < 0.1% for CMS energies ~ 200 GeV. This is summarized in 
Tables |I| and 0, in which we compare the size of the correction Ao"(i^), labeled ^4/, with 
the size of the respective NL non-ISR 0{q) correction for the 4-lepton, 2-lepton-2-quark, 
and 4-quark final states, with and without the cuts of Ref. [|lT| at 200 GeV. 

Thus, in what follows, we shall ignore Ao"(i^), as our ultimate precision tag objective. 



< 0.5%, does not require that we keep it. It will be analyzed in more detail elsewhere [10 
Further, for the cross section <y{Ya) we have already presented, for version 1.13, in Figs. 1- 
8 of Ref. |jl3|, for the csiui, i = e' , fi~ final states, the W~^'~ angular distributions in 



the e~^e~ CMS system, the W~^~ mass distributions, the distributions of the final-state 
lepton energies in the LAB frame {e~^e~ CMS frame), and the final-state lepton angular 
distributions in the W~ rest frame their corrections (relative to the Born-level). The main 
effect on these differential distributions of the improved normalization of version 1.14 is to 
shift the normalization, as we discussed above. Thus, we do not repeat their presentation 



here. We refer the reader to the results in Ref. |T3| for an investigation into the size of the 
EW=NL and FSR effects in the cases listed above insofar as YFSWW3-1.14 is concerned 
with the understanding that the shapes of the distributions apply directly to version 
1.14, but that the normalization of the EW correction should be reduced by —0.3% to 
—0.4%. In general, we found in Ref. ITBI that, depending on the experimental cuts and 



acceptances, both the FSR and the EW corrections were important in precision studies of 
these distributions; this conclusion still holds for version 1.14, of course. For example, in 
the lepton decay angle distribution, for the BARE acceptance (the final charged lepton is 
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not combined with any photons), where both the FSR and the EW correction modulate 
the distribution, whereas for the CALO acceptance of Ref. |jl3[ (all photons within 5° of 
the final-state charged lepton are combined with it) the FSR effect is almost nil whereas 
the EW correction effect remains at the level of ~ 2.0%. Here, we focus on the total CMS 
photon energy distribution (Fig. la), and the CMS photon angular distribution (Fig. lb). 
We show these results both for the BARE and CALO acceptances, as defined in the 4 



fermion Section of the Proceedings of the LEP2 MC Workshop ||TT 



In Fig. la, we see that the total photon energy distributions are different for the BARE 
and CALO cases but that the NL non-ISR correction does not affect them strongly. In 
Fig. lb, we see that, for both the BARE and CALO cases, the NL correction does affect the 
photon angular distribution away from the beam directions, as we expect. Note that this 
is the NL correction implied by the YFS exponentiation of our exact 0{a)prod correction. 
Finally, in Fig. 2, we show the effect, in the W mass and angular distributions, of using 
the screened Coulomb correction from Ref. 0], as against the usual Coulomb correction 



from Ref. [^. The effect we see is a 5 MeV shift in the peak position, associated with the 
difference between the screened and usual Coulomb corrections; we see almost no effect, as 
expected, associated with this difference on the e~^e~ CMS W angular distribution. Since 
we calculate the finite pt n{'y) corrections to these distributions, these results are new. 



Indeed, in Ref. |TT| it is shown that the results from Racoon WW and YFSWW3 (Best) 



for the distribution in Fig. la differ by ~ 20% and, as we have the dominant 0[o?^ LL 



corrections to this distribution whereas in Fig. 20 in Ref. [|Tl| the RacoonWW result only 
has the exact 0[oi) Born result for the hard photon observable, we expect that most of 
this discrepancy would be removed if the dominant O(a^) LL corrections were included in 



the RacoonWW results. This has recently been confirmed in Ref. |3^, where the authors 
of RacoonWW show that, when they include the latter corrections in their predictions 
for the cos d^ distribution in Fig. la, the discrepancy is reduced to the level of < 5%. In 



summary, from the results in Ref. [|13| and those presented here, we see that the FSR 
and EW corrections are necessary for a precision study of the distributions in the PF-pair 
production and decay process at LEP2 energies. 



We have made a detailed comparison between our results and those from Ref. ||16 
based on the program RacoonWW in the context of the LEP2 MC Workshop [|1T|. A 
complete unpublished preliminary description of the respective results of this comparison 
has appeared in Ref. ||TT|. Here, we focus on the normalization comparison of the two 
calculations at LEP2 energies. We show in Table 3 the comparison of the RacoonWW 
and YFSWW3-1.14 results for the cross sections as indicated, without cuts at 200 GeV 
(we have looked at the lower energies 184 and 189 GeV and the comparisons there are 
similar, if not better). In Table 4, we show the analogous comparisons with the LEP2 



MC Workshop cuts as described in Ref. [11 



We see that for all channels considered, the two sets of results agree to the level of 0.3%. 
This gives a total precision estimate of 0.4% for the theoretical uncertainty on the 200 GeV 
CMS energy WW signal cross-section normalization when allowance is made for further 
possible uncertainties in the higher-order radiative corrections and the implementation of 



the LPA [jTT|. This is a significant improvement over the originally quoted ~ 2% for this 
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uncertainty when the NL non-ISR 0{a) corrections are not taken into account [^. An 
effort to further reduce this 0.4% is in progress. 

Finally, with an eye toward the LC projects, we have made simulations using YFSWW3- 
1.14 for a CMS energy of 500 GeV. We show our results in Table H for the total cross section 



without cuts; here we again compare them to the corresponding ones from RacoonWW 116 
The NL corrections are significant in these results. Precision studies at LC energies must 
take these effects into account. As expected, the percentage difference between YFSWW3- 
1.14 and RacoonWW remains below 0.5% at 500 GeV CMS energy and is somewhat larger 
than at 200 GeV CMS energy. 

In summary, we have presented two recipes for combining YFSWW3 and KoralW- 
1.42 to arrive at a gauge-invariant calculation of the lyVT-pair production and decay in 
which the YFS-exponentiated exact 0{a)prod corrections are taken into account as well as 
the 0{a'^) LL FSR and YFS-exponentiated 0{a^) LL ISR correction to the background 
processes. We have illustrated our calculation with several sample MC results and we 
have compared our results on the cross-section normalization with those on Refs. [jT^] at 
200 GeV. In this way, new precision tag of 0.4% has been established for this normalization, 
which represents a considerable improvement over the original result [^ of ~ 2% when 
NL non-ISR corrections are not taken into account. 
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No cuts 


aww [fb 


(54/ %] 


5^iv [%] 


Final state 


Program 


Born 


ISR 


Born 


ISR 


V^^^T'V^ 


YFSWW3 
KoralW 


219.793(16) 
219.766(26) 


204.198(09) 
204.178(21) 


0.041 


0.044 


-1.92(4) 


(Y-K)/Y 


0.01(1)% 


0.01(1)% 


ud^^u^ 


YFSWW3 
KoralW 


659.69(5) 
659.59(8) 


635.81 (3) 
635.69 (7) 


0.073 


0.073 


-1.99(4) 


(Y-K)/Y 


0.02(1)% 


0.02(1)% 


udsc 


YFSWW3 
KoralW 


1978.37(14) 
1977.89 (25) 


1978.00 (09) 
1977.64(21) 


0.060 


0.061 


-2.06(4) 


(Y-K)/Y 


0.02(1)% 


0.02(1)% 



Table 1: The total WW cross sections aww = (^{K^), ^(Kt) at the Born and ISR level, the 4/ 
corrections ^4/ = Aa{K)/aBorniy) and the 0{a) NL correction 6^\y = Aa{Ya) / aBomiY) 
at EcM = 200 GeV. The numbers in parentheses are the statistical errors corresponding to 
the last digits of the results. All of the results are without cuts. 



With cuts 


aww [fb] 


6^f [%] 


xNL rq/1 

Oww [M 


Final state 


Program 


Born 


ISR 


Born 


ISR 


Uf,fi+T'Ur 


YFSWW3 
KoralW 


210.938(16) 
210.911(26) 


196.205 (09) 
196.174(21) 


0.041 


0.044 


-1.93(4) 


(Y-K)/Y 


0.01 (1)% 


0.02(1)% 


udfi'i'^ 


YFSWW3 
KoralW 


627.22 (5) 
627.13(8) 


605.18(3) 
605.03 (7) 


0.074 


0.074 


-2.00 (4) 


(Y-K)/Y 


0.01(1)% 


0.02(1)% 


udsc 


YFSWW3 
KoralW 


1863.60(15) 
1863.07(25) 


1865.00(09) 
1864.62(21) 


0.065 


0.064 


-2.06 (4) 


(Y-K)/Y 


0.03 (2)% 


0.02(1)% 



Table 2: The total WW cross sections aww = (^(-^s); cr(^) at the Born and ISR level, the 
4/ corrections ^4/ = Aa{K)/aBorniY) and 0{a) NL correction 5^^ = Aa{Ya) / aBomiY) 
at EcM = 200 GeV. The numbers in parentheses are the statistical errors corresponding to 
the last digits of the results. All of the results are with the bare cuts of Sect. 4.1 of Ref. [in]. 
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No cuts 


CTtot [fb] 


final state 


program 


Born 


best 


U^fl+T-U^ 


YFSWW3 
RacoonWW 


219.770(23) 
219.836(40) 


199.995(62) 
199.551(46) 


(Y-R)/Y 


-0.03(2)% 


0.22(4)% 


udjj," Ufj, 


YFSWW3 
RacoonWW 


659.64(07) 
659.51(12) 


622.71(19) 
621.06(14) 


(Y-R)/Y 


0.02(2)% 


0.27(4)% 


udsc 


YFSWW3 
RacoonWW 


1978.18(21) 
1978.53(36) 


1937.40(61) 
1932.20(44) 


(Y-R)/Y 


-0.02(2)% 


0.27(4)% 



Table 3: Total cross-sections for CC03 from RacoonWW and YFSWW3 at ^i = 200 GeV 
without cuts. The numbers in parentheses are statistical errors corresponding to the last digits. 



With hare cuts 


o"tot fb] 


final state 


program 


Born 


best 


u^^+T-Vr 


YFSWW3 
RacoonWW 


210.918(23) 
211.034(39) 


192.147(63) 
191.686(46) 


(Y-R)/Y 


-0.05(2)% 


0.24(4)% 


udjj,^ u^ 


YFSWW3 
RacoonWW 


627.18(07) 
627.22(12) 


592.68(19) 
590.94(14) 


(Y-R)/Y 


-0.01(2)% 


0.29(4)% 


udsc 


YFSWW3 
RacoonWW 


1863.40(21) 
1864.28(35) 


1826.80(62) 
1821.16(43) 


(Y-R)/Y 


-0.05(2)% 


0.31(4)% 



Table 4: Total cross-sections for CC03 from YFSWW3 and RacoonWW at v^i = 200 GeV 
with bare cuts of Sect. 4.1 in Ref. [|ll| (see the text). The numbers in parentheses are 
statistical errors corresponding to the last digits. 
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No cuts 


o-tot [fb] 


Final state 


Program 


Born 


Best 


z/^/i+r"z/^ 


YFSWW3 
RacoonWW 


87.087(11) 
87.133(23) 


89.607(32) 
90.018(27) 


(Y-R)/Y 


-0.05(3)% 


-0.46(5)% 


udfi-Uf, 


YFSWW3 
RacoonWW 


261.377(34) 
261.400(70) 


279.086(97) 
280.149(86) 


(Y-R)/Y 


-0.01(3)% 


-0.38(5)% 


udsc 


YFSWW3 
RacoonWW 


783.93(11) 
784.20(21) 


868.14(31) 
871.66(27) 


(Y-R)/Y 


-0.03(3)% 


-0.41(5)% 



Table 5: Total cross-sections for CC03 from RacoonWW and YFSWW3 at v^ = 500 GeV 
without cuts. The numbers in parentheses are statistical errors corresponding to the last digits. 



(a) 



10^ 

102 

10« 
10-1 
10-2 

10-=* 



r dil^E. 



10 



-4 



10 



-5 



da 



1) 






• ISR: BARE - ;lo4 

* Best: BARE ] 
° ISR: CALO - 



'»» 



o Best: CALO 



*"«*.. 



'»« 



*♦♦«« 



»««< 



*♦»♦ 



»». 



'»♦», 



♦♦•« 



»*«. 



IQS 



10^ 



IQi 



SE^ [GeV] - 10" 

1 

50 -1.0 



da 
d cos 6 



7 



[fb 



(b) 



• ISR: BARE 

* Best: BARE 
° ISR: CALO 
o Best: CALO 



****jj..j.. ......jjj*** 



******************* 



** *****§^?°°ooooooooooooooo°°°^*?** 






COS^^ 



10 20 30 40 50 -1.0 -0.5 0.0 0.5 1.0 

Figure 1: Distributions of the total photon energy (a) and cosine of the hardest photon polar 
angle (b) for the ud + jjiv^ + 71(7) final-state. The solid, open circle, star, and diamond 
curves correspond to the LL BARE, LL CALO, 0{a)prod 0{a^) LL FSR BARE and CALO 
YFS-exponentiated results, respectively. 
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BARE 



BARE 



U.iU 


if 


_ da-daBorn * (1): ISR + CC : 


0.05 


d(^Born (2): ISR + SCC : 
• (1) - (2) ^ 




: f 


t 

^K - 

n : 


0.0 




ft.,,. 


ffftfttttTTtTT||TTTtTtTtTTTtTTtTtTTTTTfTTT,, -^ 


0.05 




*** 






**** - 

*«* : 

*** : 


0.10 


- 


cos ^i;[,'- ; 



0.10 



0.05 



0.0 



-0.05 



-0.10 



-0.15 



: ^ _ da-daBorn : 
li, daBorn ': 

Bit 1 

; 11^ M^'-M|^'~5Mey ; 


1 III lite 


UtUit 


: 






; \ [ 




; M_[GeV] ; 


" 1 1 1 " 



-1.0 -0.5 0.0 0.5 1.0 77.5 80.0 82.5 85.0 

Figure 2: Effects of the screened Coulomb correction (SCC) on the distributions of the polar 
angle (left) and the invariant mass (right) of W^ in comparison with the usual Coulomb 
correction (CC) at Ecm = 200 GeV. As indicated the star, solid diamond and large dot 
curves are the ISR + usual Coulomb correction, ISR + screened Coulomb correction and their 
difference respectively, in the presence of YFS-exponentiation. Results are for the e+e^ — > 
W~^W~ — > ud^^u^ channel. The bare cut is that of Sect. 4.1 of Ref. pl[| . 



